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Was  the  spread  of  agropastoralism  from  the  Fertile  Crescent  throughout  Europe  influenced  by  extreme 
climate  events,  or  was  it  independent  of  climate?  We  here  generate  idealized  climate  events  using 
palaeoclimate  records.  In  a  mathematical  model  of  regional  sociocultural  development,  these  events 
disturb  the  subsistence  base  of  simulated  forager  and  farmer  societies.  We  evaluate  the  regional  simu¬ 
lated  transition  timings  and  durations  against  a  published  large  set  of  radiocarbon  dates  for  western 
Eurasia;  the  model  is  able  to  realistically  hindcast  much  of  the  inhomogeneous  space-time  evolution  of 
regional  Neolithic  transitions.  Our  study  shows  that  the  consideration  of  climate  events  improves  the 
simulation  of  typical  lags  between  cultural  complexes,  but  that  the  overall  difference  to  a  model  without 
climate  events  is  not  significant.  Climate  events  may  not  have  been  as  important  for  early  sociocultural 
dynamics  as  endogenous  factors. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Between  10  000  and  3000  cal  BC,  western  Eurasia  saw  enorm¬ 
ous  cultural,  technological,  and  sociopolitical  changes  with 
the  emergence  of  agropastoralism,  permanent  settlements,  and 
state  formation  (Barker,  2006).  Human  population  experienced 
a  dramatic  increase  (Bocquet- Appel,  2008 ;  Gignoux  et  al„  2011 ),  and 
people,  plants  and  animals  moved  or  were  moved  great  distances 
(e.g.,  Zohary  and  Hopf,  1993). 

While  the  Holocene  possibly  defines  the  start  of  major  anthro¬ 
pogenic  global  environmental  change  (Lemmen,  2010;  Kaplan  et  al., 
2011),  it  also  marks  the  period  where  climatic  shifts  could  have 
affected  human  subsistence  more  severely  than  ever  before: 
reduced  mobility  after  investments  in  settlement  infrastructure 
most  likely  increased  the  sensitivity  of  the  novel  farmers  to  envi¬ 
ronmental  alterations  (Janssen  and  Scheffer,  2004).  There  remains, 
however,  considerable  uncertainty  on  whether  and  how  climate 
instabilities  had  influenced  the  development  and  spread  of  agro¬ 
pastoralism  in  Eurasia  (Berglund,  2003;  Coombes  and  Barber,  2005). 

1.1.  Origin  and  spread  of  western  Eurasian  farming 

The  Neolithic  originated  most  probably  in  the  Fertile  Crescent, 
between  the  Levantine  coast  and  the  Zagros  ridge.  In  this  region, 
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almost  all  European  food  crops  and  animals — wheat,  barley,  cattle, 
sheep,  pigs — had  been  domesticated  and  inserted  into  a  broad 
spectrum  of  foraging  practices  during  the  tenth  millennium  cal  BC 
(Flannery,  1973;  Zeder,  2008).  Neolithic  (farming  based)  life  style 
emerged  not  before  the  9th  millennium  BC  in  this  core  region 
(Rosen  and  Rivera-Collazo,  2012),  and  expanded  to  Cyprus  by 
8500  cal  BC  (Peltenburg  et  al.,  2000);  around  7000  cal  BC,  agro¬ 
pastoralism  appeared  on  the  Balkan  and  in  Greece  (Perles,  2001). 
Propagating  in  a  generally  northwestern  direction,  agropastoralism 
finally  arrived  after  4000  cal  BC  on  the  British  isles  and  throughout 
northern  Europe  (Sheridan,  2007);  in  a  western  direction,  the 
expansion  proceeded  fast  along  the  Mediterranean  coast  to  reach 
the  Iberian  peninsula  at  5600  cal  BC  (Zapata  et  al.,  2004). 

1.2.  Transitions  and  climate 

It  has  been  argued  that  a  precondition  of  agriculture  was  the 
relatively  stable  environment  of  the  Holocene  (Feynman  and 
Ruzmaikin,  2007),  and  that  only  in  this  stable  environment  active 
cultivation  and  establishment  of  infrastructure  such  as  fields  and 
villages  was  favored  (van  der  Leeuw,  2008).  Within  its  relative 
stability,  however,  the  Holocene  climate  exhibited  variability  on 
many  spatial  and  temporal  scales  with  pronounced  multi¬ 
centennial  and  millennial  cycles  (Mayewski  et  al.,  2004;  Wanner 
et  al.,  2008).  In  addition,  non-cyclic  anomalies  have  been  identi¬ 
fied  (Wirtz  et  al.,  2010),  most  prominently  the  so-called  8.2  and  4.2 


66 


C.  Lemmen,  K.W.  Wirtz  /  Journal  of  Archaeological  Science  51  (2074)  65 —72 


events  (around  6200  and  2200  cal  BC,  respectively,  von  Grafenstein 
et  al.,  1998;  Cullen  et  al.,  2000).  Although  the  regional  scale  and 
intensity  of  the  4.2  event  has  been  strongly  questioned  (e.g.,  Finne 
et  al.,  2011 ),  the  event  had  evoked  the  formulation  of  hypotheses  on 
the  connection  between  climatic  disruptions  and  societal  collapse 
(Weiss  et  al„  1993;  DeMenocal,  2001).  Similarly,  the  globally 
documented  8.2  event  has  been  linked  to  the  abandonment  of 
many  settlements  in  the  Near  East  and  simultaneous  appearance  of 
new  village  structures  in  southeast  Europe  (Weninger  et  al.,  2005). 

It  might  be  coincidental  that  the  8.2  and  4.2  events  define  the 
time  window  of  the  Neolithic  expansion  in  Europe,  but  the  general 
view  that  environmental  pressure  on  early  Neolithic  populations 
may  have  stimulated  outmigration  has  been  put  forward  since  long 
(Childe,  1942).  Dolukhanov  (1973),  Gronenborn  (2009,  2010),  or 
Weninger  et  al.  (2009)  suggest  that  climate-induced  crises  may 
have  forced  early  farming  communities  to  fission  and  move  in  order 
to  escape  conflicts.  Berger  and  Guilaine  (2009),  to  the  contrary,  see 
the  role  of  climate  events  rather  in  creating  opportunities:  the  rapid 
farming  expansion  into  the  Balkan  could  have  been  stimulated  by 
an  increase  of  natural  fires  after  the  8.2  event,  which  opened  up  the 
formerly  forested  landscape. 

1.3.  How  sensitive  was  the  Neolithization  to  climate? 

The  relevance  of  climate  variability  and  external  triggers  for 
prehistoric  agricultural  dynamics  has  been  severely  questioned 
(e.g.,  Erickson,  1999;  Coombes  and  Barber,  2005).  Alternative 
theories  of  the  Neolithic  transition  underline  the  agency  of  early 
societies  (Shanks  and  Tilley,  1987;  Whittle  and  Cummings,  2007). 
On  the  other  hand,  the  development  of  technological,  social,  and 
cultural  complexes  can  hardly  be  thought  to  evolve  independently 
of  their  variable  environments;  and  the  spatio-temporal  imprint  of 
the  Neolithization  in  Eurasia  requires  a  geographic  approach  which 
resolves  how  people  and/or  goods  and  practices  migrated  over  long 
distances.  Berglund  (2003),  e.g.,  suggested  a  stepwise  interaction 
between  agriculture  and  climate  but  found  no  strong  links  for 
northwest  Europe. 

The  dispersal  of  agriculture  into  Europe  has  long  been  mathe¬ 
matically  formulated  based  on  Childe’s  (1925)  observation  on  the 
spatio-temporal  distribution  gradient  of  ceramics  that  Ammerman 
and  Cavalli-Sforza  (1971)  formulated  as  the  ‘wave  of  advance’ 
model.  This  simple — and  also  the  later  more  advanced  ones 
(Ackland  et  al.,  2007;  Galeta  et  al.,  2011;  Davison  et  al.,  2006) — 
diffusion  models  received  support  from  linguistic  (e.g.  Renfrew, 
1987)  and  archaeogenetic  work  (e.g.  Balaresque  et  al.,  2010).  The 
dispersal  of  agriculture  in  these  models  occurs  concentrically,  and 
can  be  modulated  by  topography  and  geography.  This  dispersal 
model  is  not  able  to  describe  the  inhomogeneous  spatio-temporal 
distribution  of  radiocarbon  dates,  which  are,  e.g.,  apparent  in 
regionally  different  stagnation  periods  (‘hypothese  arythmique’, 
Guilaine,  2003;  Rasse,  2008;  Schier,  2009). 

Stagnations  are  visible  in  the  simulation  by  Lemmen  et  al.  (2011 ), 
who  integrate  endogenous  regional  sociocultural  dynamics  with 
the  dispersal  of  agriculture.  Their  approach  connects  social 
dynamics — as  optimally  evolving  agents — to  regionally  and 
temporally  changing  environments;  in  addition,  they  account  for  the 
spatio-temporal  spread  of  populations  and  technological  traits.  Their 
Global  Land  Use  and  technological  Evolution  Simulator  (GLUES)  has 
proven  to  produce  realistic  hindcasts  of  the  origin  and  distribution  of 
agropastoralism  and  concomitant  cultures  around  the  globe  (Wirtz 
and  Lemmen,  2003;  de  Vries  et  al.,  2002),  for  Eastern  North  Amer¬ 
ica  (Lemmen,  Unpublished),  the  Indus  valley  (Lemmen  and  Khan,  in 
press),  and  western  Eurasia  (Lemmen  et  al.,  2011 ).  Using  GLUES  and 
a  globally  synchronous  climate  forcing  signal,  Wirtz  and  Lemmen 
(2003)  found  a  general  delay  of  the  simulated  regional  Neolithic 


due  to  climate  fluctuations;  at  a  global  scale,  differences  in  hind- 
casted  sociocultural  trajectories  proved  to  be  largely  independent  of 
temporal  disruptions. 

We  here  use  temporal  disruptions  that  are  defined  as  excursions 
of  a  climate  variable  far  from  the  local  mean  climate,  i.e.  extreme 
climate  events;  we  do  not  consider  rapid  climate  shifts  that 
abruptly  alter  the  climate  mean  state  (e.g.,  Dakos  and  Scheffer, 
2008).  The  hypothesis  that  extreme  climate  events  had  significant 
impacts  on  the  Neolithization  of  Europe  is  critically  examined:  we 
employ  GLUES  as  a  deductive  tool  to  reconstruct  the  Neolithic 
transition  in  Europe  and  evaluate  the  simulated  reconstruction 
against  the  radiocarbon  record  of  Neolithic  sites  in  two  experi¬ 
ments:  (1)  one  including  climate  events,  represented  by  a  pseudo- 
realistic  spatially  resolved  climate  event  history  for  the  period 
9500—3000  cal  BC;  and  (2)  another  without  climate  events. 

2.  Material  and  methods 

2.1.  Reconstructing  climate  event  history 

We  used  a  data  collection  of  134  globally  distributed,  high- 
resolution  (<200  a)  and  long-term  (>4000  a)  palaeoclimate  time 
series  collected  from  public  archives  and  published  literature.  The 
collection  only  contains  studies  where  the  respective  authors 
indicated  a  direct  relation  to  climate  variables  such  a  precipitation, 
temperature,  or  wind  regime  (e.g.  Bond,  1997;  Wick  et  al.,  2003; 
Chapman  and  Shackleton,  2000;  Gasse,  2000).  A  large  part  of  this 
data  set  (122  time  series)  was  previously  analyzed  by  Wirtz  et  al. 
(2010)  for  extreme  events;  a  complete  overview  of  time  series  in 
this  collection  is  provided  in  the  supporting  online  material 
(Table  SI).  Due  to  the  different  types  of  proxies  originating  from 
both  marine  and  terrestrial  sites  (mostly  5lsO,  see  Table  1)  the 
relation  to  climate  variables  is  often  ambiguous,  also  in  sign.  This 
ambiguity  does  not  affect  our  analysis,  as  we  are  only  interested  in 
the  spatio-temporal  characterization  of  extreme  events:  a  drastic 
excursion  from  a  climate  mean  state  stressed  regional  habitats  and 
human  populations  regardless  of  its  direction. 

Our  data  set  comprises  134  palaeoclimate  time  series,  all  of  them 
long-term  and  high-resolution,  and  provides  the  best  spatial  and 
temporal  coverage  of  any  study  we  are  aware  of.  Previous  collec¬ 
tions  used  18,  50,  60  or  80  records  (Wanner  et  al.,  2008;  Mayewski 
et  al.,  2004;  Holmgren  et  al.,  2003;  Finne  et  al.,  2011,  respectively), 
mostly  limited  to  the  last  6000  years.  The  coverage  we  use  here  is 
sufficient  to  represent  climate  variability  in  almost  all  land  areas  of 
the  world  (with  sparsest  regional  coverage  in  central  Australia, 
Saharan  Africa,  and  Northern-Central  Eurasia),  considering  the 
spatial  coherence  of  climate  signals  within  1500  km  distance  found 
by  Wirtz  et  al.  (2010)  for  their  similar  data  set. 

From  the  global  data  set,  26  time  series  are  located  in  or  near  our 
focus  area  western  Eurasia  (Table  1).  For  these  time  series  we 
analyzed  the  non-cyclic  event  frequency  according  to  the  proce¬ 
dure  in  Wirtz  et  al.  (2010):  time  series  were  detrended  with 
a  moving  window  of  2000  years  and  smoothed  with  a  moving 
window  of  50  years,  then  normalized  (Fig.  lb).  Events  were 
detected  whenever  a  time  series  signal  exceeded  a  confidence 
interval  with  threshold  p  =  1  —  1/n,  where  n  is  the  number  of  data 
points  (Thomson,  1990),  and  where  each  event  is  preceded  or  fol¬ 
lowed  by  a  sign  change  in  the  time  series. 

For  each  simulation  region,  events  from  spatially  overlapping  or 
nearby  proxy  locations  were  used  to  construct  an  aggregated  event 
time  series  specific  to  this  region  (Fig.  la— c):  (1)  a  Gaussian  filter 
with  t  =  175a  (corresponding  to  the  dating  uncertainty  in  many 
records)  standard  deviation  was  applied  to  each  event;  (2)  the 
distance  of  the  proxy  location  to  the  simulation  region  was  used  to 
assign  exponentially  decreasing  weights  to  each  event;  (3)  all  event 
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Table  1 

Table  of  26  palaeoclimate  time  series  and  associated  timing  of  extreme  climate 
events.  This  is  a  sub-set  (relevant  for  the  Western  Eurasian  focus  area)  of  the  entire 
global  collection  used  to  generated  extreme  climate  events;  the  full  data  set  is 
shown  in  the  supplementary  material  Table  SI.  Abbreviations:  Lk  =  Lake,  Cv  =  Cave, 
SST  =  sea  surface  temperature,  T  =  temperature,  T7  =  July  temperature, 
P  =  precipitation,  GSD  =  grayscale  density,  Ti  =  Titanium  content,  and 
HSG  =  hematite  stained  glass. 


No 

Site 

Proxy 

Events  (cal  BC) 

Reference 

49 

Trop  Atlantic 

SST 

None 

DeMenocal 
et  al.,  2000 

50 

N  Atlantic 

HSG 

4100 

Bond  et  al.,  2001 

52 

Ireland 

5lgO 

8960/8140 

McDermott,  2004 

53 

NW  Morocco 

5lgO 

7850/3930/4080 

Hodell  et  at,  2001 

54 

SW  Europe 

AT 

8970/3320 

Davis  et  al.,  2003 

55 

NW  Europe 

AT 

None 

Davis  et  al.,  2003 

57 

Swiss  Alpes 

T7 

7950/6030 

Wick  et  at,  2003 

58 

Swiss  Alpes 

P 

7920 

Wick  et  at,  2003 

59 

Swiss  Alpes 

P 

7920/5980 

Wick  et  at,  2003 

60 

C  Italy  Mg/Ca 

T 

5950/3270 

Drysdale  et  al.,  2006 

62 

C  Italy 

o180 

5950 

Drysdale  et  al.,  2006 

63 

S  Germany 

S,80 

8130/8050/7190 

von  Grafenstein 
et  al.,  1998 

67 

Sweden 

GSD 

9180/8350 

7700/3760 

Rubensdotter  and 
Rosqvist,  2003 

68 

N  Finland 

T 

6940/5160 

Husum  and 

Hald,  2004 

69 

N  Atlantic 

GSD 

9380/9300 

8090/3790 

Chapman  and 
Shackleton,  2000 

70 

Marmara  Sea 

S,80 

3590/3230 

Sperling  et  al.,  2003 

74 

Black  Sea 

Ti 

7220 

Bahr  et  al.,  2005 

75 

SE  Europe 

AT 

8370/7820 

Davis  et  al.,  2003 

76 

Kola 

5lsO 

8180/3550 

Jones  et  al.,  2004 

81 

Soreq  Cv, 

Israel 

51sO 

7170/4450 

Bar-Matthews 
et  al.,  1999 

82 

Red  Sea 

S,80 

9090/3730 

Seeberg-Elverfeldt 
et  al.,  2004 

93 

NW  Siberia 

T7 

None 

Hantemirov  and 
Shiyatov,  2002 

129 

Lk  loannina, 
Greece 

5lsO 

4050 

Frogley  and 

Griffiths,  2001 

132 

Lk  Acigol, 
Turkey 

S,80 

3160 

Roberts  et  al.,  2008 

133 

Lk  Frassino, 
Italy 

o,80 

None 

Roberts  et  al.,  2008 

134 

Lk  Golhisar, 
Turkey 

5lsO 

3070 

Roberts  et  al.,  201 1 

densities  were  summed  and  filtered  with  a  running  mean  with 
window  t;  (4)  all  events  representing  single  anomalies  above  the 
mean  event  density  were  used  for  further  analysis,  irrespective  of 
their  magnitude.  On  average,  seven  events  were  detected  per  time 
series. 

From  this  analysis,  we  obtain  a  pseudo-realistic  database  of 
spatially  and  temporally  resolved  climate  events.  By  using  this 
idealized  approach,  we  show  a  way  to  overcome  issues  raised  by, 
e.g.,  Schulting  (2010)  on  chronological  resolution  and  spatial 
representation  problems  of  individual  records.  In  our  study,  we 
consider  the  impact  of  the  number  of  extreme  events  and  their 
spatial  patterning  rather  than  single  chronologies.  This 
idealized — or  potential — climate  events  database  thus  allows  us  to 
go  forward  with  analyzing  the  climate— human  relationship  while 
the  reliability  of  individual  palaeoclimate  reconstructions  is  still 
being  questioned. 

2.2.  Global  Land  Use  and  technological  Evolution  Simulator 

The  Global  Land  Use  and  technological  Evolution  Simulator 
(GLUES,  Wirtz  and  Lemmen,  2003;  Lemmen,  2010;  Lemmen  et  al„ 
2011)  was  developed  to  study  how  differences  in  cultural  trajec¬ 
tories  at  the  region  scale  can  be  attributed  to  the  specific  adaptation 
of  local  societies.  GLUES  mathematically  resolves  the  dynamics  of 


Fig.  1.  Generation  of  idealized  climate  events,  (a,  top)  Map  with  example  region  and 
nearby  locations  (colored  triangles)  where  long-term  and  high-resolution  palae¬ 
oclimate  information  is  available;  concentric  circles  indicate  the  distance  from  the 
simulation  region,  (b)  Associated  palaeoclimate  time  series,  detrended  and  normal¬ 
ized;  for  each  time  series,  the  confidence  level  p  is  indicated  and  visualized  as  the 
width  of  the  background  shading.  Climate  extremes  are  identified  outside  this  confi¬ 
dence  interval  and  are  highlighted,  (c,  bottom)  Contribution  of  individual  extreme 
events  from  palaeoclimate  time  series  (weighted  by  distance  and  color-coded  as 
above)  to  the  probability  distribution  of  event  occurrence  in  the  simulation  region.  The 
dotted  line  shows  the  number  of  time  series  contributing  to  the  event  generation  for 
this  region. 

local  human  populations'  density  and  characteristic  sociocultural 
traits  in  the  context  of  a  changing  biogeographical  environment. 
One  of  the  characteristic  traits  represents  available  technologies, 
a  second  one  the  share  of  farming  and  herding  activities,  and  the 
third  one  the  number  of  established  agropastoral  economies.  Trait 
characteristics  adapt  according  to  a  growth-benefit  gradient 
dynamics  (e.g.  Smith  et  al„  2011,  for  ecological  applications).  Traits 
are  further  exchanged  between  simulation  regions  by  information 
dispersal  or  migration.  The  model  is  described  in  more  detail  in  the 
supporting  online  material. 

Our  climate  event  data  base  is  used  to  impose  sudden  reduc¬ 
tions  (f)  of  the  land  utility  whenever  a  climate  event  occurs  in 
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a  simulation  region.  If  a  local  simulated  population  experiences 
reduced  food  yield  during  a  climate  event,  people  attempt  to 
mitigate  the  crisis  by  moving  to  adjacent  regions.  In  this  respect, 
climate  events  act  as  a  stimulus  to  migration  in  the  model.  If  out¬ 
migration  is  not  possible,  population  density  declines;  with  severe 
population  decline,  also  some  of  the  technologies  are  lost  (Lemmen 
and  Wirtz,  2010). 

The  simulation  is  started  at  9500  sim  BC.1  All  of  the  685  bio- 
geographically  defined  regions  are  initially  set  with  farming  activity 
at  4%  and  established  agropastoral  communities  at  25%,  what 
represents  a  low  density  Mesolithic  technology  population  and 
a  broad  spectrum  foraging  life  style  with  low  unintentional  farming 
activity.  We  define  a  local  population  as  Neolithic  when  the  share  of 
agropastoralists  is  larger  than  the  share  of  foragers — regardless  of  its 
technology,  economic  diversity,  or  population  density  (Fig.  2). 
Experiments  are  performed  with  different  impact  strength  (f)  of 
extreme  events  on  land  utility  to  assess  climate  related  sensitivities; 
here  we  compare  the  simulation  without  climate  events  (J=  0%)  to 
simulations  where  the  utility  reduction  was/=  10%— 100%;  results 
are  shown  for/=  40%,  a  value  that  represents  a  moderate  but  not 
excessive  impact. 


2.3.  Reference  data 

Our  reference  is  a  sub-set  of  the  comprehensive  data  collection 
by  Pinhasi  et  al.  (2005),  who  used  site  data  provided  by  the  United 
Kingdom  Archaeology  Data  Service,  the  Central  Anatolian  Neolithic 
e-Workshop  (CANeW),  the  radiocarbon  CONTEXT  database,  and 
the  Radiokarbondaten  Online  (RADON)  database.  In  their  compi¬ 
lation,  Pinhasi  et  al.  included  only  sites  with  small  dating  uncer¬ 
tainty  (<200  a)  and  report  dates  based  on  calibration  of  original 
14C  measurements  with  CalPal  2004.  The  data  set  contains  dates 
distributed  across  Western  Eurasia  and  was  recommended  by  its 
authors  for  comparison  with  simulation  studies  (see  Fig.  3  for  dates 
and  locations).  For  numerical  comparison  with  the  simulation,  we 
calculate  the  timing  of  the  Neolithic  onset  from  the  radiocarbon 
date  of  all  sites  located  within  each  simulation  region.  As  there  is  no 
standard  procedure  for  calculating  such  an  area-averaged  onset 
(Muller,  2000),  we  us  the  average  date. 

3.  Results 

In  the  reference  GLUES  simulation  including  climate  variability, 
farming  originates  in  the  Fertile  Crescent  and  southeast  Europe  in 
the  7th  millennium  sim  BC  and  penetrates  into  Europe  in  a  north¬ 
west  direction  (Fig.  3).  By  3500  sim  BC  all  of  continental  Europe  has 
converted  to  farming  as  the  predominant  subsistence  style. 

3.1.  Expansion  of  agropastoralism 

The  initial  development  progresses  slowly  and  at  a  low  level.  It 
begins  during  the  67th  century  sim  BC  in  the  Levant  and  Greece,2 
followed  by  the  central  Balkan  (66th  century).  From  these 
centers,  farming  spreads  to  the  western  Black  Sea  coast  during  the 
62nd  century,  and  is  present  throughout  the  entire  Balkan,  Ana¬ 
tolia,  and  Mesopotamia  by  the  59th  century.  In  these  areas,  the 
northward  spread  of  farming  stagnates,  while  farming  activity 
intensifies  up  to  the  56th  century  on  the  Balkan,  and  while  farming 
expands  to  Italy  and  the  eastern  Black  Sea  coast. 


1  we  use  the  age  scale  'simulated  time  BC'  (sim  BC)  to  distinguish  between 
empirically  determined  age  models  (cal  BC)  and  the  simulation  time  scale. 

2  We  use  current  geographic  names  to  refer  to  the  simulation  regions. 
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Fig.  2.  Simulated  climate  events  and  evolution  of  farmer  share  for  an  example  region, 
(a,  top)  Declining  resource  availability  at  times  where  the  palaeoclimate  record  indi¬ 
cates  climate  extremes,  (b,  bottom)  Transition  from  hunting-gathering  to  agro¬ 
pastoralism  in  the  simulation  without  (gray)  and  with  (black)  climate  events.  Both  the 
timing  of  the  transition  as  well  as  the  duration  of  the  transition  (defined  as  the  time 
needed  in  the  simulation  to  increase  the  farmer  share  from  4.1%  to  95%)  are  prolonged 
by  the  introduction  of  climate  events  in  this  example  region. 

From  the  55th  to  the  59th  century  sim  BC,  a  rapid  expansion 
through  central  Europe — or  what  is  archaeologically  seen  as  the 
Linearbandkeramik  area — is  simulated,  such  that  farming  is  the 
major  subsistence  style  throughout  central  and  southeastern 
Europe,  including  the  northern  Black  Sea  coast.  From  the  46th 
century,  farming  emerges  along  the  Baltic  Sea  coast,  and  has 
become  the  dominant  subsistence  style  throughout  the  North 
European  plains,  Poland,  and  Denmark  by  the  42nd  century.  In 
France  and  England,  first  farming  is  evident  from  the  39th  century. 

Independently,  farming  also  originates  in  North  Africa  at  the 
Gibraltar  Strait  in  the  61st  century  sim  BC,  from  where  it  spreads 
into  Morocco  and  the  Iberian  peninsula,  penetrating  southern 
Spain  and  Morocco  by  the  54th  century.  The  expansion  into 
northern  Iberia  stagnates  until  the  49th  century;  and  slowly 
connects  along  the  Mediterranean  coast  to  the  Fertile  Crescent 
expansion  branch  by  the  46th  century.  All  of  the  Iberian  peninsula 
has  converted  to  farming  by  the  42nd  century.  By  the  34th  century, 
all  of  continental  western  Eurasia  and  England  rely  almost  exclu¬ 
sively  on  farming.  An  animation  showing  the  expansion  of  agri¬ 
culture  into  Europe  is  shown  as  a  Supplementary  Movie  SI. 

Supplementary  video  related  to  this  article  can  be  found  at 
http://dx.doi.org/10.1016/jjas.2012.10.023. 

While  the  simulation  captures  the  eastern  route  of  the  Neolithic 
into  Europe  (via  Greece,  the  Balkan,  Hungary,  then  north-  and 
westward,  see  Rasse,  2008  for  an  overview  of  routes)  quite  well,  it 
fails  to  simulate  the  Mediterranean  route  (Cyprus,  Sicily,  French 
coast,  then  northward);  this  was  attributed  by  Lemmen  et  al.  (2011 ) 
to  the  lack  of  sea  transport  in  the  current  model.  A  western  route 
into  Europe  is  suggested  by  the  model  emanating  from  the  Strait  of 
Gibraltar;  apart  from  some  zooarchaological  evidence  (Anderung 
et  al.,  2005),  this  route  has  not  been  confirmed  by  archaeology 
(Gronenborn,  2009;  Rasse,  2008). 

The  regional  timing  of  agropastoralism  is  contrasted  with  the 
median  radiocarbon  dates  of  Neolithic  sites  compiled  by  Pinhasi 
et  al.  (2005)  (Fig.  3).  From  this  synoptic,  time-integrated  perspec¬ 
tive,  the  simulated  centers  of  agropastoralism  in  the  Fertile  Cres¬ 
cent,  in  northern  Greece  and  at  the  Strait  of  Gibraltar  are  evident,  as 
well  as  the  southeast  to  northwest  temporal  gradient  of  the 
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Fig.  3.  Timing  of  the  transition  to  agropastoralism  in  Western  Eurasia.  The  simulated  transition  (background  pastel  shading)  is  contrasted  with  the  radiocarbon  ages  of  Neolithic 
sites  from  Pinhasi  et  al.  (2005,  solid  color  triangles).  Bold  lines  indicate  regional  stagnation  periods  where  the  onset  lag  between  neighboring  regions  is  at  least  500  years. 


Neolithic  transition.  The  general  pattern  of  the  simulated  transition 
resembles  the  pattern  that  can  be  seen  from  the  radiocarbon  dates; 
locally,  many  dates  deviate  from  the  large-scale  simulation  or 
disagree  with  each  other  within  a  simulation  region. 

Overall,  the  model  skill  visualized  in  Fig.  4  indicates  that 
recorded  variability  both  between  and  within  European  regions  can 
be  sufficiently  well  reproduced  by  the  combination  of  migration 
and  endogenous  dynamics  as  formulated  in  GLUES.  There  is 
a  significant  correlation  between  the  hindcasted  onsets  and  the 
reconstructed  onset  in  both  the  simulation  with  and  without 
climate  events  (r2  =  .37,  r2  =  .43,  n  =  39,  p  >  .99,  respectively).  The 
simulation  with  climate  events  hindcasts  the  onset  for  most  regions 
162  years  later  than  the  average  reconstructed  onset;  the  mean 
model  bias  is  =400  ±  1000  years;  the  difference  to  the  simulation 
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Fig.  4.  Reconstructed  timing  of  agropastoralism  versus  simulated  onset  using 
a  constant  or  a  fluctuating  climate  forcing  (gray  and  black  circles,  respectively).  The 
lower  left  cluster  of  outliers  represents  Fertile  Crescent  and  Anatolian  founder  regions. 


without  climate  events  (median  onset  282  years  earlier  than  the 
reconstruction)  is  statistically  not  significant. 

On  average,  climate  events  delay  the  onset  by  461  ±  176  years. 
Small  delays  (minimum  65  years)  occur  in  regions  where  climate 
events  are  temporally  separated  from  the  transition,  while  long 
delays  (maximum  1150  years)  occur  in  regions  where  several 
climate  events  occur  shortly  before  and  in  the  initial  phase  of  the 
transition.  Climate  events  during  the  transition  tend  to  prolong  the 
duration  of  the  transition  by  50—100  years  but  this  trend  is  not 
statistically  significant. 

4.  Discussion 

4.3.  Stagnation  lines  in  model  and  archaeological  evidence 

The  Global  Land  Use  and  technological  Evolution  Simulator  is 
able  to  hindcast  a  realistic  spatio-temporal  pattern  of  the  intro¬ 
duction  of  farming  and  herding  into  Europe  between  7000  and 
3500  sim  BC.  Simulated  transitions  toward  agropastoralism 
compare  well  to  a  large  data  set  of  radiocarbon  dated  Neolithic 
sites.  In  the  simulation,  as  well  as  in  the  data,  agropastoralism  did 
not  expand  uniformly,  but  rather  in  periods  of  rapid  spread  inter¬ 
rupted  by  periods  of  spatial  stagnation — but  local  intensification. 
This  rapid  Neolithization,  for  example,  applies  to  the  expansion 
from  Greece  to  the  central  Balkan  in  the  67th  century  sim  BC,  which 
is  followed  by  a  400  a  period  of  relative  stagnation.  A  very  similar 
pattern  is  hindcasted  for  the  Linearbandkeramik-like  Neo¬ 
lithization  in  the  55th  and  54th  century  cal  BC,  and  for  the  relative 
stagnation  before  the  onset  of  a  Funnelbeaker-like  cultural  complex 
further  north,  in  the  49th  to  46th  century.  Other  stagnation  lines 
identified  in  the  model  occur  around  the  Levant  and  along  the 
Pyrenees,  and  divide  the  southern  from  the  northern  Iberian 
peninsula. 

These  stagnation  periods  were  archaeologically  recognized  in 
Guilaine’s  (2003)  ‘hypothese  arythmique’.  Based  on  this  hypothesis, 
Rasse  (2008)  identified  stagnation  lines  by  analyzing  isolines  in  the 
European  Neolithic  radiocarbon  record:  he  also  found  long  dura¬ 
tion  stagnation  lines  running  through  Iberia  and  the  Pyrenees,  and 
one  crossing  Anatolia;  shorter  stagnation  lines  were  identified  in 
Greece,  one  bisecting  the  Balkan,  another  running  along  the  Car¬ 
pathian  mountains,  and  one  along  the  Alps.  Extending  the  earlier 
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work,  Schier  (2009)  emphasized  a  major  stagnation  line  along  the 
northern  edge  of  the  European  loess  belt,  separating  the  Line- 
arbandkeramik  and  Funnelbeaker  cultures  between  5100  and 
4400  cal  BC. 

Dispersal  models  of  the  European  Neolithic  create  spatial 
structure  in  the  simulated  Neolithization  pattern  from  geographic 
and  topographic  constraints  (e.g.  Davison  et  al„  2006;  Ackland 
et  al.,  2007;  Galeta  et  ah,  2011).  We  add  to  this  infrequent  climate 
excursions  emerging  as  sudden  decreases  of  the  natural  produc¬ 
tivity  and  diffusion  rates  of  traits  and  people.  This  diffusion  not  only 
depends  on  the  topography  and  geography,  but  again  on  (simu¬ 
lated)  people  and  their  technology.  Compared  to  those  dispersal 
studies,  our  model  draws  a  more  heterogeneous  picture  of  the 
regional  transitions  to  agriculture;  it  is  the  only  model  which  is 
capable  of  representing  stagnations  in  the  expansion.  We  are  thus 
one  step  further  in  establishing  a  simulation  of  the  deterministic 
pattern  of  the  Neolithic  in  Europe.  Where  this  deterministic  pattern 
fails  to  reproduce  archaeological  information,  we  can  in  the  future 
more  precisely  identify  times  and  regions  where  human  agency 
played  a  dominant  role. 

Before  the  discrepancies  can  be  exploited  scientifically,  several 
caveats  about  the  simulation  have  to  be  considered:  (1)  the  model 
deviation  is  larger  than  the  uncertainty  associated  with  radio¬ 
carbon  dates  for  most  sites.  (2)  each  simulation  region  covers 
a  large  area  which  cannot  be  fully  represented  by  individual  sites, 
such  that  the  scale  difference  introduces  an  additional  comparison 
uncertainty;  (3)  there  is  a  spatial  bias  in  the  data:  the  number  of 
sites  per  simulation  region  varies,  and  there  is  good  coverage  along 
the  transect  from  the  Levant  to  northwestern  Europe,  but  few  or  no 
information  on  eastern  Europe,  on  the  Iberian  peninsula,  and  in 
North  Africa. 


4.2.  Role  of  climate  events 

Reconstructed  climate  events  episodically  depressed  population 
density  and  to  a  lesser  extent  technological  stages.  The  impact  of 
events  on  the  timing  of  the  Neolithic  simulation  is  largest  when¬ 
ever  they  occurred  before  or  during  the  first  stage  (between  10  and 
50%  agropastoralists)  of  a  local  transition:  In  29  regions  (out  of  the 
53  simulation  regions  in  Europe),  a  climate  event  occurred  during 
the  transition.  On  average,  the  transition  was  delayed  by  50—100 
years  (w.r.t.  a  simulation  without  events).  In  three  cases,  a  region¬ 
ally  emerging  agropastoral  life  style  reverted  back  to  hunting  and 
gathering  during  a  climate  event. 

Not  everywhere,  however,  where  a  climate  event  coincides  with 
the  transition,  there  is  a  delay;  and  vice  versa,  there  are  delays  in 
regions  which  did  not  experience  climate  events  during  the  tran¬ 
sition  to  agriculture,  such  as  in  most  parts  of  western  central 
Europe.  All  individual  regions  are  closely  connected  by  neighbor 
and  remote  trade  relations:  impacts,  such  as  those  of  climate 
events,  may  have  had  far-reaching  consequences  beyond  the  locally 
affected  region  (Weninger  et  al.,  2005). 

As  evident  from  comparing  this  study  with  the  one  by  Lemmen 
et  al.  (2011),  climate  fluctuations  hindered  the  development  in 
Greece  and  the  Balkan  (which  is  now  more  realistic),  while  the 
temporal  gaps  between  the  central  Balkan,  the  Linearbandkeramik 
culture,  and  the  Funnelbeaker  regions  become  more  pronounced, 
which  is  also  more  consistent  with  archaeological  evidence 
(Pinhasi  et  al.,  2005;  Barker,  2006;  Schier,  2009).  Compared  to  the 
field  evidence,  the  simulation  shown  here  is  too  late  in  reproducing 
the  transition  to  farming  in  France.  Though  migration  is  not 
enhanced  by  climatic  triggers  in  continental  Europe,  it  does 
increase  so  for  Great  Britain,  where  the  transition  occurs  as  early  as 
France  thanks  to  increased  immigration. 


The  differences  we  obtain  between  simulated  timing  and  the 
radiocarbon  age  of  related  sites  (Fig.  3)  indicate  a  partial 
improvement  over  the  simulation  without  climate  events  shown  by 
Lemmen  et  al.  (2011):  the  overall  bias  of  the  simulation  increased 
slightly  (the  event  forcing  led  to  an  overall  delay  of  the  onset),  but 
the  variability  in  the  data  could  be  redrawn  more  precisely.  Despite 
the  non-significant  regressions  between  reconstructed  and  hind- 
casted  onsets  and  durations  of  the  Neolithic,  the  scatter  alignment 
is  close  to  the  1:1  line— apart  from  a  bias  in  the  onset  hindcasts  for 
the  Levante  regions  (Fig.  4);  neither  regionally  clumped  chronolo¬ 
gies  nor  model  assumptions  can  be  expected  to  be  quantitatively 
precise. 

Overall,  uncertainties  in  the  simulated  timing,  the  radiocarbon 
dates,  and  in  the  regional  up-scaling  are  yet  so  large,  that  the 
differences  between  the  simulation  with  and  without  climate 
events  are  not  statistically  significant  for  Europe.  Yet,  the  recon¬ 
struction  of  climate  events  improved  the  hindcasted  Neolithic 
transition,  both  in  timing  and  regional  transition  durations.  On  the 
one  hand,  this  points  to  an  effective  albeit  limited  sensitivity  of 
societal  dynamics  to  environmental  shifts.  On  the  other  hand,  this 
result  supports  the  notion  of  a  large  resilience  of  cultures,  just  like 
Grosman  (2003)  found  for  the  Late  Natufians.  Agency,  the  capacity 
of  cultural  complexes  to  adapt  and  preserve  their  life-styles  still 
appears  as  the  dominant  model  for  understanding  past  human 
ecodynamics. 

This  short  paper  is  only  a  start  in  the  direction  of  tackling  the 
climate— culture  link  on  large  scales,  and  it  gives  a  perspective  of 
what  could  be  achieved  by  bringing  palaeoclimate  data,  cultural 
modeling,  and  chronologies  together.  The  idealization  of  climate 
events  we  present  here  needs  to  be  discussed  and  evaluated 
further,  likewise  the  (few)  assumptions  made  by  the  model.  In  this 
attempt,  the  increasing  availability  of  local  chronologies  will 
stimulate  regional-scale  simulation  studies  for  Europe  and  other 
continents. 

5.  Conclusion 

We  presented  a  spatially  explicit  mathematical  model  of  the 
Neolithization  of  western  Eurasia  from  7000  BC  to  3500  BC.  Our 
model  incorporates  endogenous  socio-technological  dynamics,  as 
represented  by  the  adaptation  of  characteristic  population  traits 
and  their  interaction  with  demographics  and  changing  environ¬ 
ments.  The  simulation  improved  in  hindcasting  the  Neolithic 
expansion  after  integrating  quasi-realistic  climate  events,  espe¬ 
cially  with  regard  to  stagnation  phases  visible  also  in  field  data.  This 
outcome  helps  to  confine  the  relevance  of  climate  variability  in 
influencing  past  human  ecodynamics  and  suggests  a  dominance  of 
endogenous  over  exogenous  control  factors. 
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